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Abstract 

The e + e~ — > 7r 7r°7 process was studied in the SND experiment 
at VEPP-2M e + e~ collider in the energy region 0.60-0.97 GeV. From 
the analysis of the energy dependence of measured cross section the 
branching ratios B(u — > 7r°7r°7) = (6.6lg g ^ 0-6) ■ 10~ 5 and B(p — > 
7r vr 7) = (4llJ;g ± 0.3) • 10~ 5 were obtained. 
PACS: 13.65. +i,14.40.Cs 



1 Introduction 

In 1998 and 2000 the experiments with Spherical Neutral Detector (SND)fl] 
at VEPP-2M e + e~ collider were carried out in the energy range E = 360 — 
970 MeV where cross section of e + e~ annihilation into hadrons is determined 
by the p and u> meson decays. The integrated luminosity of 9 pb _1 collected 
in the experiment corresponds to 3.5 -10 6 and 7T0 6 produced p and oj mesons, 
respectively. One of the goals of the experiment was the investigation of the 
rare process 

e + e~ — >• p, uj — > 7r°7r°7. (1) 
* e-mail: druzhinin@inp.nsk.su 
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Figure 1: Two mechanism of p — > 7r°7r°7 decay: (a) is the VMD mechanism, (b) 
is the transition through pion loop. 

Our preliminary study of the process (|l|) was based on 1/3 of collected 
statistics. Its results were the first measurement of B(p — ► 7r 7r°7) = (4.8lf;|± 
0.2) • 10~ 5 and the measurement of B(u -> 7r°7r°7) = (7.8 ± 2.7 ± 2.0) ■ 10~ 5 
confirming the only previous measurement of this decay by GAMS: B(u — ► 
vr vr 7) = (7.2 ±2.5) • 10~ 5 §. 

The theoretical study of the p, uj — > 7r 7r°7 decays was begun by P.Singer 
in Ref. ||] where the transitions via u;7r (Fig. |I|a) and p°7r° intermediate states 
were suggested. The vector meson dominance (VMD) calculation with these 
intermediate states leads to branching ratios ~ 1 • 10~ 5 and ~ 3 • 10 -5 for 
p — > 7r°7r°7 and uj — > 7r 7r°7, respectively 0. 

For p — > 7r°7r°7 decay another mechanism through the pions loops (Fig. |TJd) 
is also possible || . The branching ratios expected for this mechanism in dif- 
ferent models || [7], |], ^ are listed in Table 0. It was noted in Ref. [§] that 
the p — ► 7r°7r°7 decay via chiral loops can be interpreted as p —>■ (77, where 
cr is a scalar state decaying into tttc pair. The dependence of B(p — > 7r°7r°7) 
on a parameters was studied in Ref. 0. The range of B(p — > 7r°7r°7) values 
in LaM model corresponds to different o widths. Since the amplitudes of 

Table 1: The branching ratios for different contributions to p — * 7r°7r°7 decay. The 
chiral loop or p — > (77 contribution is calculated in the framework of Chiral Pertur- 
bation Theory {\PT) (6|, [5J , Unitarized Chiral Perturbation Theory (U X-PT) 0] , 
Linear Sigma Model (LaM) Ji] and a pole model . 



model 


VMD 


loops, (77 


total 


XPT 
U X PT 
LaM 
a pole 


1.3 ■ 10~ 5 
1.5 • 10~ 5 
1.3 • 10~ 5 
1.0 • 10" 5 


1.0 ■ 10~ 5 
1.5 ■ 10" 5 
(0.8-2.1) ■ 10~ 5 
~ 2 ■ 10" 3 


2.9 ■ 10~ 5 
4.2 ■ 10~ 5 
(2.8-4.7) ■ lO" 5 
~ 2 • 10" 3 
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p — ► dJ7r° — > 7r°7r°7 transition (Fig. |l|a) and of the pion loops (Fig. |l|b) are 
of the same order of magnitude, their interference is substantial. The in- 
terference contribution into branching ratio is predicted to be positive. The 
theoretical values for the total branching ratios are also listed in Table |l|. 
The predictions of the chiral models || [7|, §, |J are in agreement with our 
previous experimental result. The significantly larger value of B(p — > 7r°7r°7) 
was obtained in Ref. [|l(| using a pole model. Their result contradicts to 
existing experimental data. 

In u — > 7r°7r°7 decay the contribution of pion loops is G-parity suppressed 
while the contribution of kaon loop is small due to large kaon mass. There- 
fore, it is assumed that the u — > 7r°7r°7 decay proceeds through p°7r° inter- 
mediate state. The first measurement B(u — ► 7r 7r°7) = (7.2 ± 2.5) • || 
significantly exceeded the existent prediction of VMD model: 3 ■ 10~ 5 ||. An 
attempt to explain this discrepancy was done in Ref. [IT] where p-u mixing 
was taken into account and coupling constants were extracted from exper- 
imental values of T{u — > 37r) and r(p° — > 7r°7). As a result the estimated 
value of B(u — ■> 7r°7r°7) increased up to (4.6 ± 1.1) • 10~ 5 . Similar results 
(4.5-4.7)-10~ 5 were then obtained in Refs. || ||. In Ref. [0 the large ex- 



perimental value of B{u — > 7r°7r°7) was explained by additional contribution 
of the uj — > o"7 transition and used to extract the value of g wul coupling 
constant. 

In the present work we present the experimental results on B(u — ► 7i 7c 0/ ~f) 
and B(p — > 7r 7r°7) based on full SND data sample. 



2 Event selection 

For analysis five-photon events with the energy deposition in the calorimeter 

E tot > 0.7 • E (2) 

and the total momentum measured by the calorimeter 

P tot < 0.15 -E/c (3) 

were selected. Here E is e + e~ center of mass energy. 

Due to high beam background rate in 5% of events fake photons ap- 
pear. This makes possible for lower photon multiplicity QED processes 
e + e~ — > 27, 37,47, and p, uj — > n°j,ru — > 37 decays to imitate five-photon 
events producing main background contribution for the process under study. 
Detector response to the beam background was studied using special events 
recorded with a random generator trigger. The information on the fired de- 
tector channels in these events was used for simulation of the process under 
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study and the background processes. Considerable suppression (by a fac- 
tor of 8) of the background from events with fake photons was achieved by 
imposing the following cuts: 

E min > 30 MeV, 30° < min < 150°, (4) 

where E min and 6 min are the energy and polar angle of the softest photon 
in an event. These cuts reduce the detection efficiency for the process under 
study by 25%. Another background source is the e + e~ — > 777 — > 37r°7 — ► 7-y 
reaction producing five-photon events mainly due to the merging of near pho- 
tons. To suppress this background, the parameter x-y describing transverse 
energy deposition profile of the detected photon [O was used. The cut 



X 7 < 5 (5) 

suppresses the e + e~ — > 777 background by a factor of 2 with a 5% loss of 
actual 5-photon events. 

Further selection was based on the kinematic fitting of the events. Com- 
patibility of the event kinematics with e + e~ —>■ 57 and e + e~ — ► 37 hypotheses 
was checked. For the 37 hypothesis two out of five photons were considered 
spurious: all 3-7 subsets were tested and the best one with minimum \ 2 
value was selected. As a result of kinematic fitting the \ 2 values, X57 an d 
X3 7 , were calculated for both hypotheses. The cut 

Xs 7 > 20 (6) 

practically eliminates e + e~ — > 27, 37 background with the loss only 2.5% of 
the events of the process under study. Figure |2| depicts the X57 distribution 
of the experimental and simulated events. The following cut was imposed on 
this parameter: 

Xs 7 < 20. (7) 

Finally, the events with two 7r° mesons were selected. To do this the kine- 
matic fit in e + e~ — > 7r°7r°7 hypothesis was performed and the following cut 
was imposed: 

X7T7T7 X5 7 < 10. (8) 

Here Xwiry is the \ 2 value of the kinematic fit for the e + e~ — > 7r 7r°7 hypoth- 
esis. The Xmr-y — X5-y distributions for the experimental and simulated events 
are shown in Fig. [3]. 

The difference between distributions of the background events and events 
of the process under study (Figs. || and |3|) was used to estimate the accuracy 
of the background calculation. The experimental distribution was fitted by 
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Figure 2: The %5 7 distribution. Points 
with error bars represent experimental 
data. The histogram is a simulation of 
the process (jlj) and background pro- 
cesses. The shaded histogram shows 
the contribution of background pro- 
cesses. 



Figure 3: The Xnwy~ X57 distribution. 
The points with error bars are exper- 
imental data. The histogram is the 
sum of simulated events of the process 
(||) and background processes. The 
shaded histogram is the contribution 
of the background processes. 



the sum of simulated distributions for the process (|l|) and background. As a 
result of the fit the ratio K = Nj£g/N%%. = 0.7 ± 0.2 was determined. For 
softer selection criteria without cuts this ratio is equal to 1.30 ± 0.05 

and does not depend on the beam energy at our level of statistical accuracy. 
From this we conclude that the systematic error of the background estimation 
does not exceed 50%. 

The total of 310 events were selected with estimated background of 15 ±7 
events. The main sources of the residual background are e + e~ — > 777 and 
e + e~ — > 47 processes. The distribution of selected events and calculated 
background over center of mass energy is given in the Table 0. No events 
were detected below 600 MeV. The uncertainty of the center of mass energy, 
integrated luminosity, and detection efficiency are listed in Table || for each 
energy point. The uncertainty of the center of mass energy includes the 
beam energy spread and the energy shift between 1998 and 2000 scans. The 
integrated luminosity was measured using e + e~ — ► 77 process. The statistical 
error of the luminosity in each energy point does not exceed 1% and is not 
included in the table. Its systematic error was estimated to be 3%. The 
detection efficiency for the process (0) was determined by simulation. The 
differential cross section of the e + e~ — > 7r°7r°7 process calculated in VMD 



model WM was used for simulation. The systematic error of the detection 
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efficiency, including the model error due to possible contribution from o~7 
intermediate state was estimated to be 5%. 



3 Fitting of the cross section 

The fitting procedure maximizes the logarithmic likelihood function 

L = £lnP(iVf p ,Af), 



where Pj is a Poisson probability to detect observed number of events N^ xp 
in the i-th energy bin with a theoretical expectation of Nf h . The theoretical 
expectations were calculated as 

Af =e l L l a(E l )(l + 5(E l )) + N^ 9 , 

where N^ kg is a calculated number of background events, £j is a detection 
efficiency, Li is the integrated luminosity, u(E) is e + e~ — > 7r°7r°7 cross section 
depending on a set of approximation parameters, 5 is a radiative correction. 
The radiative correction, which is a functional of the cross-section energy 
dependence cr(E) was determined within the fitting procedure. The 

values of radiative correction evaluated for each experimental energy point 
are listed in Table |2|. The model error of the (1 + 8) value does not exceed 
3%. The values of the experimental cross section calculated as 

fJ ex P _ ]\f bk 9 
exp Jv i Jv j 



e i L i {l + S{E i )) 



are shown in Fig.^ and listed in Table 0. The systematic error of the cross 
section is determined by the the errors of the detector efficiency, integrated 
luminosity, and radiative correction. It was estimated to be 7%. 

To calculate the cross section a(E) the amplitude of the e + e~ — > 7r°7r°7 
process was parametrized as 

= A^BWp + a x BW p , + a 2 BW p „) + f3A pal BW p + 1 A U) BW U) , (9) 







The first term in Eq.@ is the amplitude of the e + e~ p,p',p" — > ujtx 
where p' and p" are excitations of the p(770). Second and third terms are 
e + e~ — ► p(770) —>■ 0^ — > 7r°7r°7 and e + e _ — > uj — > 7r°7r°7 amplitudes. Each 
amplitude is written in a factorized form. The functions A puJn o, A pcrl , A^ 
depending on the momenta of final particles describe the dynamics of vector 
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Table 2: Center of mass energy E, its standard deviation <7g, integrated luminosity 
L, number of selected events N exp , calculated number of background events Nbkg, 
detection efficiency e, radiative correction 1 + 5 and the e + e~ — > 7r°7r°7 cross 
section a exp . 



E,MeV 


OE,MeV 


L, nb" 1 


N 

1 ' exp 


Nbkg 


e 


1 + 5 


a exp ,nb 


600.1 


0.29 


88.3 





0.1 


0.273 


0.912 


-O.OOSlHol 


630.1 


0.30 


116.1 





0.1 


0.269 


0.906 


-0.0041°;^ 


660.2 


0.25 


271.6 


2 


0.3 


0.273 


0.900 


0.0251HS 


690.2 


0.29 


167.2 


2 


0.2 


0.263 


0.895 


o.o46lffi 


720.3 


0.26 


588.5 


1 


0.8 


0.251 


0.892 


o.oo2±8:S? 


750.2 


0.32 


219.0 


3 


0.1 


0.259 


0.897 


0.0571°;^ 


760.2 


0.31 


238.9 


2 


0.3 


0.251 


0.896 


o.032l°:°£ 


764.2 


0.32 


250.4 


5 


0.2 


0.254 


0.892 


o.o85±8:8S 


770.2 


0.31 


284.4 


8 


0.3 


0.253 


0.877 


n 1 oo+0.063 


774.2 


0.34 


217.1 


7 


0.2 


0.252 


0.855 


A -| 4 r+0.081 


778.1 


0.34 


247.9 


6 


0.5 


0.261 


0.820 


a 1 a 4 +0.068 
U.1U4_ 045 


780.2 


0.35 


319.5 


16 


1.2 


0.263 


0.807 


a 9-1 O+0.075 


781.1 


0.33 


339.6 


20 


1.2 


0.267 


0.807 


n 9 r 7 +0.076 
U.ZO 1 -0.061 


782.1 


0.31 


656.3 


34 


1.0 


0.257 


0.815 


0.24018:82 


783.2 


0.30 


473.4 


30 


2.0 


0.253 


0.833 


0.28018:81 


784.2 


0.32 


346.2 


24 


0.7 


0.261 


0.857 


a or,-. +0.077 
U.OUl_oo63 


785.3 


0.24 


212.3 


12 


0.4 


0.257 


0.890 


a 9OO+0.094 
U.zoo_q.070 


786.1 


0.33 


267.7 


11 


0.4 


0.255 


0.914 


0.17018:8^ 


790.1 


0.34 


191.4 


4 


0.3 


0.258 


1.006 


0.07418:81^ 


794.2 


0.34 


206.7 


1 


0.2 


0.256 


1.044 


a a i 4 +0.042 
U.U1^_ 015 


800.2 


0.32 


276.8 


10 


0.3 


0.255 


1.053 


a -I oa+0.057 
U.l<JU_o.042 


810.2 


0.34 


279.5 


3 


0.4 


0.240 


1.043 


n 37 +0 - 030 

u - uo ' -0.024 


820.1 


0.36 


315.2 


2 


0.3 


0.244 


1.035 


a a 91 +0.033 
U.UZ1_ oi6 


840.2 


0.35 


677.5 


8 


0.8 


0.247 


1.025 


0.04218:8?^ 


880.0 


0.41 


376.0 


7 


0.5 


0.222 


1.001 


o.o78l8:8M 


919.9 


0.44 


478.6 


8 


0.3 


0.256 


0.916 


0.06918:81 


939.9 


0.43 


469.0 


22 


0.7 


0.248 


0.856 


a 9-, 4 +0.058 
U.Z14_ 047 


949.7 


0.32 


261.7 


20 


0.3 


0.261 


0.855 


a oqo+0.095 
U.OOO_ 076 


957.7 


0.32 


233.9 


131 


' 0.2 


0.263 


0.858 


a 949+0.089 
U.Z4/_ 067 


969.7 


0.34 


251.5 


29 


0.5 


0.250 


0.865 


a p-94+0.119 
u - Jz ^-0.099 




650 700 750 800 850 900 950 

E, MeV 

Figure 4: The cross section of the e + e~ — > 7r°7r°7 process. Points with error bars 
are experimental data. Solid line is a result of fitting in the model 3 of Table ||. 
The dashed line corresponds to the fit with B(p — > 177) = 0. 



mesons decays, 
shapes: 



The functions BWi describe the Breit-Wigner resonance 



BWi 



P, p',p",w- 



mf — E 2 — iETi(E) ' 

Here m« and Ti(E) are resonance mass and energy dependent width. The 
cross section is calculated from Eq. (Q) by integration over the phase space of 
final particles: cr(E) = J |A 7r7r7 | 2 <in. At the energy above utt threshold the 
Breit-Wigner functions of p mesons are modified BW Pi — > BW Pi C PiUJ7T , where 
CpiUiT are Blatt-Weisskopf factors, restricting fast growth of the T PiUJ7T partial 



widths [|I~6f : 



a 



a 



+ (R qiv (m Pi )y 
\ \ + {Rq u {E)) 2 



Pi 



P',p"- 



(10) 



Here is the u meson momentum in pi — > ujtt decay. The range parameter 
R is supposed to be the same for p, p', p" mesons. The main decay modes of 
p mesons were taken into account for calculation of the energy dependence of 
the resonance widths. For instance, in the case of p(770) we use the following 
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are experimental data. Curve is a fit result in the model 3 of Table pi 



expression: 

Here is a pion momentum in the p — > 27r decay. The Blatt-Weisskopf 
factor C p7r7r is expressed by the formula similar to Eg. (|TUD . 

The amplitude of the e + e~ — > ujtt process is described by formulas from 
Ref. [H, |T7[ and depends on 8 parameters. The data from energy region 



below 1 GeV are insufficient to determine them. Therefore we used additional 
measurements of the e + e~ — > ujtt cross section in the 1-1.4 GeV energy 
range by SND |I7| and of the ujtt spectral function in r — ► ujttv t decay 
by CLEO |18]. The spectral function can be converted to corresponding 
production cross section in e + e~ collisions using CVC hypothesis [jH^ . The 
all-data fit on e + e~ — > 7r°7r°7 cross section is shown in Fig.|5|. For E > 1 GeV 
experimental data are well described by three models with parameters listed 
in Table |3|. The ranges of parameter values correspond to the variation of the 
parameter R from to 2 GeV -1 . The statistical errors are not shown because 
they are significantly smaller than model biases. For models with two excited 
p states the p' mass and width were fixed to 1400 MeV and 500 MeV. These 
values are close to p' parameters from tt^tt spectral function data [2(| |21 . 
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Table 3: Parameters of the models describing e + e — > ujtt — > 7r°7r°7 cross section 
for E > 1 GeV. 







ny 




«i 






a 2 


1 


14.3-15.8 








1630-1710 


630-1000 


-(0.19-0.24) 


2 


14.1-15.7 


1400 


500 


-(0.04-0.06) 


1580-1620 


420-580 


-(0.14-0.18) 


3 


15.4-16.6 


1400 


500 


-(0.39-0.42) 


1560-1640 


380-780 


0.24-0.30 



Below ujtt threshold the amplitude of the e + e~ — > ujtt process drops 
rapidly and the product l-SW^I 2 / |A pW7r | 2 <in in contrast with the correspond- 
ing product for p — > o"7 transition does not demonstrate resonance behavior. 
This allows to separate contributions of the two p decay mechanisms by 
measurement of energy dependence of the e + e~ — ► 7r°7r°7 cross section. The 
p — ► o"7 decay amplitude was described by the xPT and LaM models from 
Ref. 0. The three sets of a parameters [^2], |23], 0] used in the LaM model 
are listed in Table [5[ The xPT model corresponds to m a — > oo. The f3 
parameter represents the difference between observed value of the p — > 0-7 
decay amplitude and theoretical prediction. 

For u; — > 7r°7r°7 decay the variation of the final state phase in the CMS 
energy interval of the u meson is small, so we cannot separate different decay 
mechanisms studying the cross section energy dependence. Thus the ampli- 
tude of the uj — > 7r°7r°7 decay was written according to VMD model ||14|| . The 
p-u mixing was taken into account following Ref. 0. Possible contributions 
of other mechanisms would result in a deviation of the complex parameter 7 
from 1. 

Full description of the energy dependence of the cross section below 1 GeV 
requires extra four parameters e.g. absolute values and phases of (3 and 7: 
|/3|, 0/3, I7I, cj)y. But we prefer two other sets of parameters: B(p — ► 0-7), <frp, 
B(u — > 7r°7r°7), 7 or B(p — > 7r°7r°7), 0/3, B(u — > 7r°7r°7), 7 . The branching 
ratios are related to (3 and 7 as 

B{p -> 7r°vr 7) = t4— / \A pLU7T (m p ) + pA pal (m p )\ 2 dU, 
B(p - <ri) = -f — / \{3A pay (m p )\ 2 cm : 



B (u -> tt°A) = 7^ — / l7^K)| 2 dn, 

where oy = 127r5(l / — > eV)/My is a total vector meson production cross 
section in e + e" collisions. 
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Table 4: The branching ratios of p and u decays (B x 10 5 ) and P(x 2 ) values 
obtained as a result of cross section fitting with different values of (f)g, <j>~. The um 
amplitude was described by the Model 3 from Table || Model 1 from Table |5] was 
used for description of p — > aj amplitude. Only statistical errors of parameters 
are shown. 



Model 


^(CJ -> 7T 7T 7) 


5(p -> 7T 7T 7) 




^(x 2 ) 


1 


0/3 » 0, 


> 7 


« 


6.3±1;| 


4 1+ 1 - 

4 - i -0.9 


1 Q+0-9 


35% 


2 


0/3 ~ 7T, </ 


i 7 


« 


12.3±?;1 


Q Q+0.9 


4.4 ± 1.0 


30% 


3 


0/3 ~ 0, </ 


) 7 


7T 


25.5^2.3 


5 1+1-0 
d - i -0.9 


1.9 ± 1.0 


6% 


4 


0/3 ~ 7T, 


>0 


« 7T 


15.8 ±2.3 


4 7+ - 9 

4 - ' -0.8 


5.6±& 


8% 



Characteristic feature of the process under study is a large interference 
between the contributions of p and u decays. For instance the cross section 
of e + e~ — > uj — > 7T 7r°7 process at = evaluated using the table value 
of B{uo -> 7r°7r°7) = (7.2 ± 2.5) • 10" 5 is equal to 0.12 nb. The interference 
with p meson increases this value up to approximately 0.25 nb (Fig.|j). The 
experimental data on the energy dependence of the cross section are insuffi- 
cient for determination of unambiguous solution for interference phases 0g, 
7 . There are four solutions listed in Table [|. The third and fourth ones 
correspond to a large destructive contribution into to decay from mechanisms 
other than p°7r°. The B(u — > 7r°7r°7) values obtained in this case disagree 
with existing experimental value B{uo — > 7r°7r°7) = (7.2 ±2.5) • 10~ 5 . The so- 
lution with (f)p m 7r, 7 ~ can be ruled out for two reasons: B(u — > 7r°7r°7) 
exceeds the table value by 1.7 standard deviations and consistency of the 
calculated spectrum of the recoil photon with the experimental one is poor. 
The analysis of the photon spectrum is described in the next section. 

For the only survivor solution with both phases close to zero, the model 
dependence of the fit parameters was studied. Three models of excited p 
states (Table |3|) and four sets of a parameters (Table |^) were tested. The 0g 
was found ranging within (20 -j- 80)° ± 80°. These values are in agreement 
with theoretically expected zero value J7|. Therefore the final fitting was 
performed with <pp = 0. The phase 7 was considered as a floating parameter 
to take into account its possible shift due to the contribution of mechanisms 
other than uj — ► pn . The fitted 7 = —(2 4- 20)° ± 20° is consistent with 
zero. 

The probabilities of the p and uj decays into 7r°7r°7 obtained with differ- 
ent parameters of a meson are listed in Table |5|. The spreads in parameter 
values correspond to different models describing e + e~ — > um° cross section 
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Table 5: The probabilities of the p and oj decays into 7r°7r°7 (B x 10 5 ) obtained with 
different parameters of a meson. The spreads in the parameter values correspond 
to the models listed in Table |3[ Theoretical values of branching ratios are taken 
from Ref. The final results with statistical and systematic errors are given in 
the bottom line of the table. 



Model 


m a 


rv 


B(u -> 7r°7r°7) 


B(p^ 


7r 7r°7) 




aj) 




MeV 


MeV 




exp. 


th. 


exp. 


th. 


LaM 


478 


324 


6.2-6.5 


4.0-4.1 


3.8 


1.5-1.9 


1.5 


LaM 


555 


540 


6.3-6.6 


4.2-4.3 


2.8 


1.9-2.3 


0.8 


LaM 


478 


263 


6.2-6.4 


4.2-4.3 


4.7 


1.7-2.1 


2.1 








6.5-6.9 


3.9-4.0 


2.9 


1.8-2.2 


1.0 




6.6±l;f ± 0.6 


4 - i -0.9 


±0.3 


1.9lg;|±0.4 



above 1 GeV. All models reproduce the experimental data well. Therefore 
parameter midrange was taken as a final result. Its spread was regarded 
as the model error. The branching ratios obtained this way with statistical 
and systematic errors are listed in the last row of Table |5|. The systematic 
error includes the model error, uncertainties in the detection efficiency and 
integrated luminosity. The variation of the background level within its sys- 
tematic error practically does not change the B(u — > 7r°7r°7) central value 
and results in following additional uncertainties of the p meson branching 
ratios: 7% for B(p — > 7r°7r°7) and 12% for B(p — ► 0-7). Since these un- 
certainties affect statistical significance of the results they were added to 
statistical errors. The energy dependence of the cross section in the model 
with rricr = 478 MeV and T a = 324 MeV is shown in Fig. |] together with 
the curve corresponding to B(p — > aj) = 0. The P(x 2 ) value for the latter 
model is equal to 0.5%. 

4 The energy and angular spectra 

From Table |2| it is seen that selected events are mainly concentrated in two 
energy regions: 180 events near u peak and 92 events in the range 920- 
970 MeV above the e + e" — ■> um° reaction threshold. The angular and energy 
distributions in the latter region agree with uir° mechanism. Recoil photon 
spectrum for events from 760-800 MeV energy range is shown in Fig.^. Al- 
though this energy region is dominated by u peak the contributions of both u 
and p decays must be taken into account to obtain the theoretical spectrum. 
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Figure 6: The recoil photon spectrum 
for the experimental events of the reac- 
tion (|]) in the energy range 760 < E < 
800 MeV (points with error bars) and 
the results of simulation in the model 
3 from Table [| (solid line) and in the 
model 2 from Table |I] (dashed line) . 



Figure 7: Distribution of cosine angle 
between photon and 7r° meson in the 
7r°7r° rest frame. Point with error bars 
are experimental data, Solid line is a 
simulation with uj decay into pure p°ir° 
state. Dashed line is a simulation with 
u) decaying into a mixture of /?°7r° and 
o"7 states. 



The spectrum calculated in model 3 from Table [| (Fig.^j) is in a good agree- 
ment with the experimental one. In this model it is supposed that uj — > 7r°7r°7 
decay proceeds through p°7r° intermediate state. Another way is to describe 
the uj decay by a sum of contributions of uj — > p°7r° and uj — > mechanisms. 
To do this we fix B(uj —>■ p°7i° — > 7r°7r°7) at 2.5T0~ 5 and fit the uj — > decay 
contribution to a value yielding observed uj — > 7r°7r°7 branching ratio. In case 
of constructive interference this leads to B(uj — > — > 7r°7r°7) = 1.3 • 10 -5 
and the photon spectrum close to that expected for p°7r° mechanism. As was 
shown in the previous section, assumption of destructive interference results 
in B(uj —* 7r°7r°7) inconsistent with the PDG table value. 

The second theoretical spectrum in Fig|] corresponds to the model 2 from 
Table |] with destructive interference of ujtt and o"7 amplitudes in p decay. 
For this model the consistency between theoretical and experimental spectra 
calculated using Kolmogorov test |3[ is about 1%, which was one of the 
reasons to discard this model. 

Additional information about mechanism of uj decay can be obtained from 
the analysis of angular distributions. One of such distribution is shown in 
Fig.[7[ The same figure displays the theoretical distributions obtained under 
assumptions that uj decay proceeds through either pure p°7r° intermediate 
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state or a mixture of p°7r° and cry. One can see that our limited statistics 
does not allow to distinguish these two models. 



5 Summary 

The branching ratios measured in this work, 

B(lo -> 7r°A) = (6.6±i;| ± 0.6) • 10~ 5 , (12) 

B(p -> 7T 7r 7 ) = (4.1±5;g ± 0.3) ■ 10~ 5 . (13) 

are in a good agreement with our preliminary results @ and GAMS mea- 
surement B(p —>■ 7rV 7) = (7.2 ± 2.5) -10 5 0, but have higher accuracy. 

The probability of p —>■ 7r°7r°7 decay significantly exceeds VMD model 
prediction (1.3 — 1.5) ■ 10~ 5 . This excess can be explained by the contribution 
of the decay via scalar state p — > cry. Two mechanisms, p — > a^y and p — > un, 
can be separated using difference in energy dependence of their amplitudes. 
Our result on p — > 0-7 decay 

B(p -> cr 7 -> 7T°A) = (1-918:1 ± 0-4) ■ 10- 5 (14) 

differs from zero by 2.4 standard deviations and is consistent with the pre- 
dictions of chiral models |7|, The magnitude of B(p — > 0-7) is sensitive to 
a parameters. As can be seen from Table [5], the models with r CT ~ 300 MeV 
give the most consistent description of the experimental data. 

The value of the branching ratio B(u — > 7r°7r°7) = (6.7 ± 1.2) ■ 10" 5 , 
obtained by averaging of our measurement with the GAMS result exceeds 
theoretical predictions, (4.6 ± 1.1) • 10~ 5 [TTJ] and (4.7 ± 0.9) • 10" 5 @, 
by 1.3 standard deviations. It is necessary to make some remarks about 



these predictions. The result of Ref. ]TTJ] is based on table value of T(p 



7r°7) = 102 ± 26 keV P5fl . It must be corrected taking into account newer 
measurement r(p — > 7i 0/ ~f) = 76 ± 22 keV [26|, which is close to the value for 



charged p, T(p ± — ► n ±r y) = 68 ± 8 keV This decreases the predicted 

B(uj — > 7r 7r°7) and worsens agreement with the experiment. In the Refs. 0] 
and H the values of g puJ7T equal to 15 GeV -1 and 15.9 GeV -1 were used to 
calculate B(u — > 7r 7r°7) oc 5 , ^ J7r fl , p7r7 oc g^ n /gf,- On the other hand, the 
use of these g pujn values for VMD calculation of T(u — > 3n) and T(lj — >• 7r°7) 
leads to too large values conflicting with experimental data. For example, the 
g P ujn obtained from uj — ► 37r width assuming intermediate p7r state is equal 
to (14.3 ± 0.2) GeV- 1 0. Therefore, our opinion is that 4.6 - 10~ 5 is the 
maximum branching ratio acceptable within VMD model and the additional 
theoretical study is required to explain the large value of B(oj — > 7r°7r°7). 
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